THE sl 3 JONES POLYNOMIAL OF THE TREFOIL: A CASE STUDY OF 

g-HOLONOMIC SEQUENCES 



STAVROS GAROUFALIDIS AND CHRISTOPH KOUTSCHAN 

Abstract. The 5(3 colored Jones polynomial of the trefoil knot is a g-holonomic sequence of two variables 
with natural origin, namely quantum topology. The paper presents an explicit set of generators for the 
annihilator ideal of this q-holonomic sequence as a case study. On the one hand, our results are new and 
useful to quantum topology: this is the first example of a rank 2 Lie algebra computation concerning 
the colored Jones polynomial of a knot. On the other hand, this work illustrates the applicability and 
computational power of the employed computer algebra methods. 
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1. The colored Jones polynomial: a case study of (7-holonomic sequences 

1.1. Introduction. The aim of this paper is to investigate the 5(3 colored Jones polynomial of the trefoil 
knot which is a g-holonomic sequence in two variables. Such sequences arise naturally in quantum topology 
when studying knots. The relevance of our results within this field become evident when taking into account 
that this is the first example of a rank 2 Lie algebra computation concerning the colored Jones polynomial 
of a knot, in accordance to the 5(3 A J conjecture of [Gar04]. 
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Using computer algebra, we compute an explicit (conjectured) set of generators for this g-holonomic 
sequence as a case study. We employ the Mathematica packages developed by M. Kauers (see [Kau09b, 
Kau09a] and the second author (see [KoulOb. KoulOa, Koull]). These computations show the power of the 
underlying computer algebra methods and may serve well the theoretical and practical needs of quantum 
topology. 

Since quantum topology and computer algebra are rather disjoint subjects, we need to briefly review some 
basic concepts before we present our results. 

1.2. q-holonomic sequences. A q-holonomic sequence (/«(?)) for n £ N is a sequence (typically of rational 
functions /«.(<?) £ Q(<?) m one variable q) which satisfies a linear recursion relation: 

ad{q n , q)fn+d{q) + ■■• + a (q n ,q)f n (q) = 

for all n£N, where aj(u, v) £ Q[tt, v] for all j = 0, . . . , d. The algorithmic significance of such sequences was 
first recognized by Zcilberger in [Zci90], where also a generalization to multivariate sequences / n (q) for n = 
(m, . . . ,n r ) £ N r was given. A down-to-earth introduction of (q-) holonomic sequences and closely related 
notions is given in [Kou09]. It is well known that q-holonomic sequences enjoy several useful properties: 

(a) Addition, multiplication, and specialization preserve q-holonomicity. 

(b) Finite (multi-dimensional) sums of proper q-hypergeomctric terms are q-holonomic in the remaining 
free variables. This is the fundamental theorem of WZ theory; see [WZ92]. 

(c) The fundamental theorem is constructive, and computer-implemented; see [PR97, PR, KoulOa, 
PWZ96]. 

1.3. The colored Jones polynomial is q-holonomic. Five years ago, a natural source of q-holonomic 
sequences was discovered: Knot Theory and Quantum Topology; sec [GL05]. Let us review some basic 
concepts of knot theory and quantum topology here. For a detailed discussion, the reader may look into 
[BN, GL05, Jan96, Jon87, Tur88, Tur94]. A knot K in 3-space is a smoothly embedded circle in the 3-sphere 
S 3 considered up to isotopy. The simplest non-trivial knot is the trefoil: 




In his seminal paper [Jon.87], Jones introduced the Jones polynomial of a knot K. The Jones polynomial of a 
knot K is an element of = %, q^ 1 ] and can be generalized to a polynomial invariant Jk.v{q) £ 2[q ] 

of a framed knot K colored by a representation V of a simple Lie algebra g. When the Lie algebra is 5(2 
and V = (D 2 is the 2-dimensional irreducible representation, J K <c 2 (?) is the classical Jones polynomial. The 
definition of Jk,v{q) uses the machinery of quantum groups and may be found in [Jan96, Tur88, Tur94]. 
Without going into details, let us point out some features of this theory. Irreducible representations of a 
simple Lie algebra g are parametrized by dominant weights, and once we choose a basis of fundamental 
weights, the set of dominant weights can be identified with N r , where r is the rank of the Lie algebra. Thus, 
the g colored Jones polynomial of a knot is a function n £ N r >-> Jx,n{q) £ Z[g ±:L ]. In [GL05] the following 
result was shown. 

Theorem 1.1. [GL05, Thm.l] For every simple Lie algebra g of rank r, and every knot K, the function 
Jk,»(q) '■ N r — > Z^ 1 ] is q-holonomic. 

It will be useful to recall in brief the proof of the above theorem using principles of q-holonomic sequences 
and quantum group theory. The quantum group invariant Jkv{q) of a knot K is a local object, obtained 

by 

• choosing a planar projection of the knot, 

• assigning a tensor (the so-called i?-matrix or its inverse, which depends on the representation V of 
the simple Lie algebra g) to each positive or negative crossing, 
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• summing over all contractions of indices, and adjust for the framing. 

It follows that Jk,v(q) is a hnite multi-dimensional sum where the summand is a product of entries of the 
i?-matrix with certain framing factors. In [GL05] it was shown that if we choose a basis of the quantum 
group suitably, the entries of the i?-matrix are g-proper hypergeometric multisums. In addition, the framing 
factor is a monomial in q raised to a quadratic form in the summation variables. It follows from property 
(b) of Section 1.2 that Jk,v(q) is g-holonomic. This summarizes the proof of Theorem 1.1. 

The proof of Theorem 1.1 is algorithmic. When the Lie algebra is 5I2 and K is a, knot with a planar pro- 
jection with d crossings, there is an explicit proper q- hypergeometric term in d + 1 variables (fci, . . . , kd, n) 
such that summation in all (ki, . . . ,kd) € N d (this is a finite sum) gives JK,n(q)- The algorithm is dis- 
cussed in detail in [GL05, Sec. 3] and has been computer-implemented in the KnotAtlas via the command 
ColouredJones; see [BN]. The above presentation shows the limitation of the WZ method. Although 
Jk,u(q) is a g-holonomic sequence presented by an explicit q- hypergeometric multisum, when K is a knot 
with several (eg. 8) crossings, the computation of a recursion relation for JK,n{q) via WZ theory of [WZ92] 
becomes impractical. 

1.4. The st3-colored Jones polynomial of the trefoil. Concrete formulas for the colored Jones polyno- 
mial Jk,v(q) are hard to find in the case of higher rank Lie algebras, and for good reasons. For torus knots T, 
Jones and Rosso gave a formula for Jt,v(q) which involves a plethysm map of V, unknown in general; see 
[RJ93]. The plethysm map was computed explicitly by [GVb] for 5(3 (and in forthcoming work of [GVa] for 
all simple Lie algebras). This gives an explicit formula for the trefoil, and more generally for the T(2,b) 
torus knots, where b is an odd natural number. Let us recall this formula, which is the starting point of our 
paper, using the notation of [GVb]. Let 

(1) fb,n u m(Q) = JT(2,b),n u n 2 {q) 

denote the 5(3 quantum group invariant of the torus knot T(2, b) colored with the irreducible representation 
Vni,n 2 °f s '3 °f highest weight A = niu>i + 7J2W2 (0-framed and normalized to be 1 at the unknot), where 
Hi,H2 are non-negative natural numbers and u>i,u>2 are the fundamental weights of 5(3; [FH91, Hum78]. 

Theorem 1.2. [GVb] For all odd natural numbers b we have 

g-2b /min{ni,re 2 }m-J b 
/&,rai,n 2 (?) = I™ 1 '" 2 X/ (~ 1 ) fc ^2ni-2fc-2i,2» 2 +fc-2i^2 2 » 1 -2k~2l,2n 9 +k-2l 

aill <" 2 \ 1=0 k=0 

min{n 1 ,n 2 } n? — l ^ 

+ (- 1 ) fc ^2«i+fc-2;.2n2-2fc-2i^2 2 n 1 +fc-2i.2«2-2fc-2i 
Z=0 fc=0 

min{ni ,712 } 

- ^2 rf 2n 1 -2i,2n 2 -2i^2 2 n 1 -2i,2™ 2 -2i 
(=0 

where the quantum integer [n], the quantum dimension d niiTl2 and the twist parameter 9 nit1l2 ofV ni . 
defined by 



(2) [n] 

(3) dfl 1,712 



g 2 - q 2 
1 _ 1 

qz — q 2 

[ni + l][n 2 + l][ni + n 2 + 2] 



[2] 



(4) nun2 = q 



-(nf+n 1 n2+n'2)+n 1 +n2 



1.5. The q-Weyl algebra. Our results below are phrased in terms of the q-Wcyl algebra, and its corre- 
sponding module theory and geometry Let us recall the appropriate q-Wcyl algebra for bivariatc g-holonomic 
sequences. 
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Consider the operators Mi, Li for i = 1,2 which act on a sequence f ni ,n 2 (q) £ Q(<7) by: 
(Lif) nu n 2 (q) = fm+l,nM) {Mif) ni . n2 (q) = q ni / ni ,„ 2 (q) 

(L2.f)n u n 2 (q) = /ni,n 2 +l(g) (M 2 f) nun2 (q) = <Z" 2 /ni ,n 2 (?) 

It is easy to see that the operators g, Mi, L\, M 2 , L 2 commute except in the following cases: 

L\M\ = qMxLx L 2 M 2 = qM 2 L 2 

It follows that if R(q, Mi,M 2 ) 6 Q(g, Mi, M 2 ) is a rational function, then we have the commutation relations 

/ = {L.Riq, Mi, M 2 ) - R(q, qM 1 ,M 2 )L 1 ,L 2 R(q, M U M 2 ) - R(q, M 1 ,qM 2 )L 2 }. 
Consider the (localized) q- Weyl algebra 

(5) W q . Xoc = QfeM^Ma)^!,^}/^) 

generated by the operators Li,L 2 over the field Mi, M 2 ), modulo the 2-sided ideal / above. W qt \ oc is an 
example of an Ore extension of the field <Q(q, Mi, M 2 ); see [Ore33]. Given a bivariatc sequence f ni ,n 2 {q) G 
(Q(g), consider the annihilating ideal 

Ann/ = {Pe W q \Pf = 0}. 

Ann/ is always a left ideal in W 9 , which gives rise to the cyclic W^joc-modulc Mj = Wq,i c/Ann/. Finitely 
generated W^joc-modules M have a well-defined theory of dimension (and Hilbert polynomial), and satisfy 
the key Bernstein inequality; see [Sab93, Thm.2.1.1]. Thus, one can define g-holonomic modules. When 
AI = Wq Xo c/J for a left ideal J, the Hilbert polynomial of the left ideal J can be computed by Grobner 
bases [Buc65] and their Grobner fans which are well-defined for the case of the Ore algebra W q j oc under 
consideration; see [Sab93, Kou09]. For a computer implementation of Grobner bases in the Ore algebra 
W qAoc , see [CS98, Kou09]. 

1.6. Our results. One may try to compute generators for the ideal of recursion relations using techniques 
which arc well known in computer algebra, but which are recalled in Section 2 for sake of self-containcdncss. 
Given elements Pi for i e S of M^.ioc, let (Pi\i G S) denote the left ideal of W qj \ c that they generate. 
Furthermore, let r denote the symmetry map that exchanges n\ and n 2 , i.e., for P S Wq^oc we have 
t(P(M 1 ,M 2 ,L 1 ,L 2 )) = P(M 2 , M U L 2 , Li). 

Theorem 1.3. (a) Let Pi, P 2 , P3 be the operators of Appendix A and 

(6) J=(P 1 ,P 2 ,P 3 ). 

Then {Pi,P2,p3} is a Grobner basis for J with respect to total degree lexicographic order (L\ > L 2 ), and 
J is a zero-dimensional ideal of rank 5. 
(b) We have 

J=(Qi,0.2) 

where {Qi,Q 2 } is a Grobner basis with respect to the lexicographic order (Li > L 2 ). The five monomials 
below the staircases of {Pi,P 2 ,P^) and (Qi,Q 2 ) are given by: 



2 






£ 2 








ii 










Li 



(c) The Grobner fan GF(J) of J is a fan in R^. with rays generated by the vectors (4, 1), (2, 1), (1, 1), (1, 2), (1,4) 
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on (Li,L 2 ) coordinates. 



(d) The left ideal J is invariant under the symmetry map r. The operator Pi is itself symmetric (modulo 
sign change) 

Proof. The computer proof of statements (a)-(d) is given in the electronic supplementary material [Koull]. 

□ 

Conjecture 1.1. Let J be the left ideal defined in (6). Then 

Remark 1.2. The colored Jones polynomial satisfies the symmetry JK,v*{q) = JK,v{q) where V* is the 
dual representation of the simple Lie algebra. Since V* i n2 = V n2 , ni for the representations of SI3, it follows 
that fb,ni,n 2 (q) = fb,n 2 ,ni (<?)• Thus, the fact that J is invariant under r gives additional evidence for the 
above conjecture. 

1.7. The 5(3 AJ Conjecture for the trefoil. One of the best-known geometric invariants of holonomic 
-D-modules are their characteristic varieties, [Kas78]. In the case of holonomic D- modules over the complex 
torus (<E*) r , the characteristic variety is a Lagrangian subvariety of the cotangent bundle T*(((D*) r ), namely a 
conormal bundle, g-holonomic Z?-modules have a characteristic variety, too. In the case of the q-holonomic D- 
modules that come from knot theory, namely the g colored Jones polynomial of a knot K , their characteristic 
variety is a subvariety of <D 2r , where r is the rank of the simple Lie algebra g. In [Gar04] the first author 
formulated the AJ Conjecture which identifies the characteristic variety of the colored Jones polynomial of a 
knot with its deformation variety. The latter is the variety of G-representations of the boundary torus of the 
knot complement that extend to G-rcprcscntations of the knot complement, where G is the simple, simply 
connected Lie group with Lie algebra q. 

In down-to-earth terms, the SI3 AJ Conjecture involves the image of the 5(3 recursion ideal of a knot under 
the partially defined map: 

(7) e : W q ,i oc — ► Q(Mi, M 2 )[Li,L 2 ], e(P) = P\ q=1 
For the ideal J of Proposition 1.3, we have 

(8) e(J) = (pi,p 2 ,p 3 ) 

where pi,P2,P3 are given in the appendix. In general, it is hard to compute the variety of SL3((D) represen- 
tations of a knot complement and to compare it with the above image. For the trefoil knot, however such a 
computation is possible (see [Gar]). It turns out that e(J) coincides with the SL3((D)-deformation variety of 
the trefoil, thus confirming the AJ conjecture. 

2. The computation 

2.1. Guessing generators for the annihilator ideal. The method of guessing is nothing else but an 
ansatz with undetermined coefficients. Assume that a multivariate sequence (/ n (g)),n £ N r , of rational 
functions in q, i.e., f n (q) S Q(g), is given by some expression that allows to compute the values of (f n {q)) 
for small n, e.g., for all n inside the hypercube [0,no] r . For finding linear recurrences with polynomial 
coefficients (they are guaranteed to exist if (/ n (?)) is known to be q- holonomic), an ansatz of the following 
form is used: 

(9) ( c a , b M a L b Wg) = 

^ (a.bJeSCN 2 '' ' 
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where the structure set S is finite and the unknowns c &i b are to be determined in Q(g). The multi- index 
notation M a = Mf 1 • • ■ M" r is employed here. Linear equations for the c a ,b are obtained by substituting 
concrete integer tuples for n into (9) and by plugging in the known values of (f n (q)). If sufficiently many 
equations are generated this way, i.e., if the resulting system of equations is overdetermined, then it is to be 
expected with high probability that the solutions of this system are indeed valid recurrences for the sequence 
in question. 

In the case of a univariate sequence (/«(<?)), usually the recursion relation with minimal order is desired. 
The analogue in the multivariate case is a set of recurrence operators that generate a left ideal in W q ^\ oc 
of minimal rank. For their many nice properties it is natural to aim at a Grobner basis [Buc65] of this 
annihilating ideal. There are two strategies how to obtain such a basis of recurrences. The first strategy 
consists in finding a bunch of recurrences by guessing and then apply Buchbcrgcr's algorithm to obtain a 
Grobner basis of the ideal they generate. The more recurrences are used as input, the higher is the probability 
of ending up with the full annihilating ideal (the one with minimal rank) . Since this latter computation can 
be quite tedious, we follow a different strategy that doesn't require Buchbcrger's algorithm: the structure 
set S in the ansatz (9) is adjusted by trial and error until the linear system for the c a ,b admits a single 
solution, which (supposedly) is an element of the desired Grobner basis. The latter is achieved by searching 
systematically (in the same fashion as in the FGLM algorithm) through the monomials L b that appear in 
the ansatz. This second strategy will finally deliver a set of operators that are very likely to be elements of 
the annihilator of (/ n (<?)) (as it is the case for the first strategy), and that are very likely to form a Grobner 
basis (which is not the case for the first strategy). 

The above described method involves a lot of trials until the structure set for each element of the Grobner 
basis has been figured out exactly. In order to make this reasonably fast, modular techniques are used. First, 
the input data (/ n (<z)) for n G [0,no] r is computed only for a specific choice of q, say q = 64, and modulo 
some prime number, say p = 2147483647. While in other examples of this type it doesn't really matter which 
value is substituted for q (there may be a small finite set of non-eligible integers only), the situation here 
is slightly different and the choice q = 64 not at random: note the rational numbers \ and i that appear 
in the exponents of Theorem 1.2. Although in the end all fractional powers disappear, yielding a Laurent 
polynomial, the do appear in intermediate expressions. Thus setting q = to 6 for some m £ N simplifies our 
job significantly. All following computations (constructing and solving the linear system) are now done with 
the homomorphic images in TL V instead of Q(<7). The computational complexity is reduced drastically since 
no expression swell can happen. Once the exact shape of (9) is found this way, the final computation over 
Q(q) can be started with a refined ansatz by omitting all unknowns c a ,b that become zero in the modular 
computation. This methodology of guessing as described above is implemented in the Mathematica package 
Guess. m (see [Kau09b, Kau09a]). 

In our case study of the sis colored Jones polynomial of the trefoil knot we have r = 2 and Theorem 1.2 
allows to compute /3,ni,n 2 (9) for sufficiently many (rii,ri2) € N 2 . Trial and error gave evidence that the 
(Li, L2)-supports of the Grobner basis elements are 



As stated in Proposition 1.3, the ideal generated by these three elements has rank 5 and the monomials that 
cannot be reduced by Pi , P 2 , P3 are the following: 



We tried to guess an operator with support U and total degree 60 with respect to Mi and M 2 in its 
coefficients, but did not succeed. This indicates that we found the ideal of minimal rank, taking into account 
that Pi, P2, P3 have (Mi, M2)-total degrees 23, 28, and 27, respectively. The modular computations involved 
solving linear systems with several thousand unknowns, which is not a big deal with today's hardware. The 
details of this section are given in the Mathematica notebook GuessIdeal.CJ.2.3.nb of [Koull]. 



supp L (Pi) 
supp L (P 2 ) 
supp L (P 3 ) 




U = {LiL 2 , Ll Li,L 2 ,l}. 



THE st 3 JONES POLYNOMIAL OF THE TREFOIL: A CASE STUDY OF q-HOLONOMIC SEQUENCES 



7 



2.2. Grobner bases in the g-Weyl algebra. Once a Grobner basis with respect to some specified term 
order is found, the FGLM algorithm (sec [FGLM93]) can be employed to obtain a Grobner basis of the same 
ideal with respect to any other term order. Using the FGLM implementation for noncommutative operator 
algebras in the Mathcmatica package HolonomicFunctions .m (see [KoulOb, KoulOa]), the second basis given 
in Proposition 1.3 was computed. Fixing a certain term order destroys the symmetry with respect to n\ 
and ri2 that is inherent in the problem of the sl3-colored Jones polynomials, because either L\ > L2 or vice 
versa. Therefore the generators Pi, P2, P3 are not expected to be symmetric themselves (Pi however is). The 
symmetry of the problem is revealed when the Grobner fan is computed. The resulting universal Grobner 
basis is symmetric with respect to n\ and ni- Moreover, specific symmetric recurrences in the ideal can be 
found by calling the command FindRelation ([KoulOb, KoulOa]). For example, by specifying the support 
{1, L1L2, L\Li2, L\L<2, LfL^} we found an operator with this support and coefficients that are symmetric with 
respect to Mi and M2. This operator gives rise to a recurrence for the diagonal sequence /3,n,n(<?)- The 
details of this section are given in the Mathcmatica notebook GroebnerBasis . CJ . 2 . 3 .nb of [Koull]. 

2.3. Changing the normalization of the colored Jones polynomial. In quantum topology, one often 
uses other normalizations of the colored Jones function of a knot. Theorem 1.2 uses the normalization where 
the 0- framed unknot has colored Jones polynomial 1. If we use the TQFT normalization, where the colored 
Jones of the unknot has value d ni ,n 2 and an arbitrary fixed framing c £ the corresponding colored Jones 
function is given by 

Pb,c,ni ,n 2 (?) = d ni ^n 2 9 ni „ 2 fb,m ,n 2 (?) 

where d ni . n . 2 and ni ,n 2 are given by Equations (3) and (4) respectively. The result is still a bivariatc sequence 
of Laurent polynomials in q. Since d nii „ 2 and ni ,n 2 are proper q-hypcrgeometric (in the language of [WZ92]), 
it follows that the generators for the annihilator ideal of fb, ni ,n 2 (<?) can easily be translated to generators of 
the annihilator ideal of P, l ni,n 2 (9) and vice- versa. This is explained in detail in the Mathcmatica notebook 
GroebnerBasis. C J. 2. 3. nb of [Koull]. 

Appendix A. Generators for the recursion ideal of f3, ni ,n 2 (q) 

P 1 = -q 6 M 3 M 2 (q 3 M 1 - l)(qM 2 - l)(q 4 M 1 M 2 - 1)(<7 5 MiM 2 - l)Fi£ 2 
~q 4 M 1 M 2 {q 2 M 1 - l)(q 2 M 2 - l)(Mi - M 2 )(q 4 A4 1 M 2 - l)F 6 L x L 2 
+q 6 M 1 Ml(qM 1 - l){q 3 M 2 - l){q A MiM 2 - l)(q 5 MiM 2 - l)r(Pi)il 
+M 2 (q 2 M 1 - l)(qM 2 - l)(q 3 M!M 2 - l)(q B M?M 2 - l)F r Li 
+ (-M 1 )(qM 1 - l)(q 2 M 2 - l)(g 3 MiM 2 - l)(g 5 M!M 2 2 - l)r{F 7 )L 2 
+q\qM 1 - l)(qM 2 - l)(Mi - M 2 )(q 2 M x M 2 - l)F 2 

P 2 = -q 21 M^M 2 7 (q 4 M 2 -l)(q 5 M 1 M 2 ~l)(q 5 M 1 Mi -l)F 1 Ll 

+q 7 M?Ml(q 2 Mi - l)(q 2 M 2 - l)(q 4 MiM 2 - l)(g 6 MiMf - l)F 8 LiL 2 
+q 6 Ml(q 3 M 2 - l)(q 4 M 1 M 2 - l)F l6 L 2 2 

+q s M?Mi(q 2 M 1 - l)(qM 2 - l)(q 3 MiM 2 - l)( 9 6 MiM 2 2 - l)F 3 Li 
+ (q 2 M 2 - l)(q 3 M!M 2 - l)( g S MiMf - 1)F 14 L 2 
-q 2 {qM 2 - l)(q 2 M 1 M 2 - 1)F 9 

P :i = -q 17 M 3 M 2 (q 2 M x - l)(g 3 M 2 - l)(q 5 M 1 M 2 - l)(g 5 AfiAf| - \)F X L X L\ 
+q 5 M 1 Mi(q 2 M 1 - l)(q 2 M 2 - l)(q 4 M 1 M 2 - l){q 6 M 1 M% - \)F X2 L X L 2 
+q 6 Ml(qM 1 - l)(q 3 M 2 - l)(q 4 M 1 M 2 - 1)F 13 L 2 2 
+q 7 M 1 Mi(q 2 M 1 - l)(qM 2 - l)(q 3 MiM 2 - l)FuLi 
-(qMi - l)(q 2 M 2 - l)( g 3 MiM 2 - l)(q 5 AhM^ - 1)F 15 L 2 
+q 2 {qKh - l){qM 2 - l)(q 2 MiM 2 - l)(g 6 MiM 2 2 - i)F 4 
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Qi = q 47 M?MZ (q 6 M 2 - l)(q 7 M 1 M 2 - l)(q e MiM^ - l)(q r MiM2 - 1)F 5 L S 2 
-q 24 M?Ml(q 5 M 2 - l)(q 6 MiM 2 - l)(g 6 MiMf - 1)F 20 F 2 
+q 9 Ml(q 4 M 2 - l){q 5 M 1 M 2 - l)( g 9 MiM 2 2 - 1)F 23 L| 
+ (q :i M 2 - l)(g 4 A/iA/ 2 - l)(g 7 MiM| - 1)F 24 L 2 
+q 2 (q 2 M 2 - l){q 3 MiM 2 - l)(g 10 AfiM| - 1)F 18 L 2 
-g 6 (gA/ 2 - \){q 2 M x M 2 - l)(g 9 MiM| - l)(g 10 MiM 2 2 - l)Fi 

Q 2 = g 8 M 1 3 M|(g-l)(g 2 Mi-l)(g 2 M 1 2 +gMi + l)(gM2-l)(g 3 MiM2-l) 
x(g 6 MiM| - l)(g 7 MiAff - l)(g 8 MiMf - l)F 5 Li 
- g 37 Mi 6 M 2 10 (g 5 M 2 - l)(g 6 MiAf 2 - l)(g 6 AfiMf - 1)F 8 £ 2 
+g 18 Af 1 2 A/ 2 6 (g 4 A/2 - l)(g s MiM 2 - 1)F 19j L| 
-q e M 2 {q 3 M 2 - l){q 4 MiM 2 - l)(g 7 MiMf - l)F 2 il| 
-(g 2 M 2 - l)(g 3 A/iA/ 2 - 1)F 22 L 2 
+q 2 (qM 2 - l)(q 2 MiM 2 - l)(g 8 MiM 2 2 - l)Fi 7 

where i 7 !, . . . , F24 are irreducible polynomials in (Q[g, Mi, M2], too large to print all of them here: 



Fi = q w MfMt - q w MtM 2 - q 15 M$M% + q 44 M\M 2 - q 14 M 4 Ml + q 13 M^M^ - q X2 M\M^. - q 14 M 3 M^ 
+q 11 M 2 Mi - q W MlMl + q 10 M 4 Ml + q 9 M 4 Ml - q 9 M 3 M 2 + q s M 4 M% + 2q 8 M 3 M 2 i + q s M 2 Mi 
-qM 4 M\ - q 7 M 3 M' 2 - q 6 M 2 M 2 - q 4 M 2 M 2 - q 4 M' 2 + q 3 M 3 + q 3 M x Ml + qM 2 - Mi 

F 2 = q 19 M\Ml - q 16 M^M 2 - g 16 A/ 4 A/ 2 6 - q 13 M?M 2 - g 13 M 4 Aff + q"M 4 M? 2 + q 11 M 3 M 2 + q w M^M% 
+q 10 M 2 Ml + q 9 M 4 M 2 + g 9 M 3 M| + q 9 M 2 M 2 - 2q 8 M 3 M 2 + q 7 M\mI - q h M 2 M 2 - q 4 M 2 M 2 
-q 4 M 1 M 2 - q 3 M 3 - q 3 M 2 M 2 - q 3 M x Ml - q 3 Ml + q 2 M 2 M 2 + q 2 MiM 2 + Mi + M 2 

F 6 = q 22 M[M 7 2 - q 20 MlM 2 - g 18 A/ 7 M 2 6 - g 18 M\ M 2 + q 15 MjM 2 + g 16 ' M^M 2 + q 15 MlM 2 - 2q K M\mI 
+q 15 M 4 M 2 7 + q 14 M^Ml + q 14 M 4 Ml + g 13 M^M| + q 2 M\m\ - q 12 M?M 2 - q 11 M?M 2 2 - g 11 M^Mf 
-q 10 M 4 M 2 + q 9 M 4 Mt + q 9 M 3 M 2 - g 8 M 1 5 Af| - 3g 8 M 4 Af 2 3 - 3g 8 M 3 Af| - q s M 2 M 2 + q 7 M 4 Ml 
+q 7 M 3 M 2 + q 6 M 3 M% + q 6 M 2 Ml + q 4 M 4 + q 4 M 3 M 2 + 2g 4 M 2 Af| + g 4 MiM 2 3 + g 4 M 2 - q 3 M 3 M 2 
-g 3 AfiAf| - gAf 2 - qM\M 2 - qM 2 + M1M2 

F 7 = q 24 M\M 7 2 - q 21 M*M 2 - q 21 M^M 2 + q 19 M^M 2 - g 18 M 7 M| - q ls M^M 2 + q u M\m\ - q 16 M 4 M 2 
+q 15 MjMl + q 15 M 4 Ml + q 14 MfMl + q 14 M^M 2 + 2g 14 A/ 4 Aif - Sq^M^M^ + q 12 MfM^ 
-q 12 M 4 M 2 + q 41 M 3 M 2 - q xl MlM\ + q 19 M\mI - g 1() M 4 M 2 3 - 2g 10 M 4 M 2 + q 9 M 3 Ml + q 9 M 2 M 2 
-q*M\M 2 + q s M 4 Mi - 3g 8 M 3 M| - q 8 M 2 M 2 + q 7 M 4 Ml + q 7 M 4 + g 7 M 3 M| - q 5 M 4 + q s MfM 2 
+q 4 M?M 2 + q 4 M? 2 - g 3 A/ 3 - q 3 M X M 2 - qM 2 + Mi 



The complete list is available at [Koull]. 



Appendix B. The q = 1 limit for the recursion ideal of /a.n^na^) 
Let pi = e(Pi) denote the image of Pi when q — 1. Then, up to polynomial factors of Mi,Mi, we have: 
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pi = Mi - L 2 M\ - LxMxM-x + L 2 M\M 2 - l\m\M 2 + LiL 2 M*M 2 - M 2 + L x Ml - m\m\ + L\m\mI 

+l 2 m\mI - l x l 2 mIm 2 + Mi Mi - L\M\m^ + l x l 2 mIm\ - l%mImI + L x M\Ml - LiL 2 M\Ml 

—L\_L 2 M\M 2 + L\MxM 2 + LiM^Mi - L\m\M 2 - L 2 M*M 2 + LlM?M 2 - L 2 M?M 2 + L x L 2 mImI 

-l^mImI + l\m\mI + l 2 m\m 2 - l\m\mI 

p 2 = (-1 + L 2 )(Mi - MiM 2 + M| - LiAff M| + L 2 M\mI + LxMlMl - L 2 MiM 2 - L 2 MtM 2 + L 2 M 2 

+L 2 M?Ml - L 2 M\Mt - LlMfMl + h\M\Ml - l\m\m\) 
P3 = {-l + L 2 )(-l + M?M 2 -MiMi+LiMiMl 

-L 2 Mi - LxMlMl + L 2 M*M 2 + L\M\M 2 - L 2 MjM 2 + L\L 2 M^M\ - L x L 2 mIm 2 + L x L 2 M\mI) 

qi = (-l + L 2 ) 2 (l + L 2 Ml)(l + L 2 M'tMl)(-l + L 2 M^Ml) 
q 2 = -(-l + i 2 )(l+L2A/|)(l + L 2 M 1 3 M 2 3 )(-l + L 2 M 1 2 M 2 4 ) 
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